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Abstract 

We discuss inclusive production of open charm mesons in proton-proton scattering at the BNL 
RHIC. The calculation is performed in the framework of /Cf-factorization approach which effec¬ 
tively includes higher-order pQCD corrections. Different models of unintegrated gluon distribu¬ 
tions (UGDF) from the literature are used. We focus on UGDF models favoured by the LHC data 
and on a new up-to-date parametrizations based on the HERA collider DIS high-precision data. 
Results of the fcf-factorization approach are compared to next-to-leading order collinear predic¬ 
tions. The hadronization of heavy quarks is done by means of fragmentation function technique. 
The theoretical transverse momentum distributions of charmed mesons are compared with recent 
experimental data of the STAR collaboration at \/s = 200 and 500 GeV. Theoretical uncertainties 
related to the choice of renormalization and factorization scales as well as due to the quark mass 
are discussed. Very good description of the measured integrated cross sections and differential 
distributions is obtained for the Jung setBO CCFM UGDF. Revised charm and bottom theoretical 
cross sections corresponding to those measured recently by the STAR and PHENIX collaborations 
for semileptonic decays of D and B mesons are presented. Significant improvement in theoretical 
description of the non-photonic electrons measurements is clearly obtained with respect to the 
previous studies within the fcf-factorization. 
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I. INTRODUCTION 


Production of heavy quarks is of substantial and ongoing interest in high energy 
hadronic collisions. This statement had not changed in nearly 40 years, when the charm 
and bottom flavoured particles were discovered. The energy scale for the production of 
charm and bottom quarks is significantly higher than the typical Quantum Chromody¬ 
namics (QCD) scale, Aqcd ^ 0-2 GeV. This gives the value of strong coupling of the 
order of as ^ 0.2 — 0.3, which is small enough to apply perturbative QCD techniques. 
At the present level of knowledge and current experimental abilities, heavy flavours are 
known as one of the best testing grounds of the theory of hard QCD interactions. Theo¬ 
retical analyses of charm and bottom cross sections in proton-proton interactions provide 
an unique precision tool in this context. Due to their mass, heavy quarks are also belived 
to be a special probe of the medium created in heavy ion collisions. Since heavy quarks 
are only produced in the initial stage of the heavy ion collisions, heavy quark distribu¬ 
tions from proton-proton interactions supply a well defined initial state. Their further 
propagation through the hot and dense medium probes its interesting characteristics. 

Measurements of charm and bottom cross sections at hadron colliders can be per¬ 
formed in the so-called indirect way. This method is based on measurement of leptons 
from semileptonic decays of open charm and bottom mesons, which are often called non¬ 
photonic. The semileptonic decay modes allow for an indirect measurement of the D and 
B meson cross sections and by further extrapolation the charm and bottom quark cross 
sections. The indirect methodologies to measure non-photonic electrons/muons have 
been used since the early 1970's IIO. 

The decay of hadrons by the weak interaction can be viewed as a process of decay of 
their constituent quarks. The charm and bottom flavours are not preserved in weak inter¬ 
actions, so their weak decays are possible. Within the semileptonic decays, the allowed 
(according to electric charge conservation) transitions & —)■ c, w and c —)• s, d are involved 
by the emission of charged W boson, which further creates a lepton and the correspond¬ 
ing antineutrino W —> Zv/. The semileptonic decay widths are proportional to the square 
of the appropriate elements of the Cabbibo-Kobayashi-Maskawa (CKM) quark-mixing 
matrix ||2t], which contains information on the strength of flavour-changing weak decays. 
The decays within the same quark generation are strongly favoured over decays between 
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generations. 

Main virtues of the semileptonic modes come from the fact that they have bigger 
branching fractions than the hadronic D meson decay channels. Moreover, in the former 
the effects of strong interactions can be isolated and thus they are better accessible exper¬ 
imentally. In addition, direct lepton production through weak and electromagnetic inter¬ 
actions in hadronic collisions is suppressed relative to the strong interaction by twelve 
and four orders of magnitude, respectively ISO. 

However, electrons and muons are certainly not rare paticles because they are abun¬ 
dantly produced in light hadron decays. These rather problematic background must be 
accounted for and eliminafed through experimental techniques. If the primary sources 
of background are well understood and substracted, the remaining events can then be 
attributed to the heavy flavour signal. 

Another method for experimental investigation of charm and bottom quark produc¬ 
tion at hadron colliders is the direct procedure based on full reconstruction of all de¬ 
cay products of open charm and bottom mesons, for instance in the 
D"*" ^ or —> }/tpK~^ —?• charmels. The hadronic decay prod¬ 

ucts can be used to built invariant mass distributions, permitting direct observation of D 
or B meson as a peak in the experimental invariant mass spectrum. Open charm D and 
B mesons are characterized by rather long lifetimes, of the order of 10 and ^ 10 
seconds, respectively Charm and bottom quarks decay essentially at the collision vertex, 
while heavy flavour mesons decay from a secondary vertex offset by the boosted life¬ 
time of the paticle. In the direct approach the charm and bottom contributions can be 
well separated, which is not the case in the indirect method. In the latter case, it can be 
achieved only within the analysis of lepton-meson (e.g. e-D) correlations 1^, which are 
easily available e.g. in the fcf-factorization approach. 

The STAR and PHENIX collaborations have measured transverse momentum distri¬ 
butions of electrons coming from the semileptonic decays of heavy flavoured hadrons 
in profon-profon scattering at the RHIC energy ^/s = 200 GeV with lepton transverse 
momenta up to 10 GeV in the midrapidity region In addition, the STAR collabo¬ 

ration was able to separate the charm and bottom contributions to the spectra of heavy 
flavour electrons ITQ. The PHENIX collaboration has also measured non-photonic dilep¬ 
ton invariant mass spectrum from 0 to 8 GeV in proton-proton collisions at ^/s = 200 
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GeV iSD. On the theoretical side, the cross sections for inclusive production of the non¬ 
photonic electrons at RHIC have been studied theoretically up to the next-to-leading or¬ 
der pQCD collinear approximation within the FONLL approach in Ref. |l9(]. The firsf the¬ 
oretical investigation within the competitive QCD fcf-factorization framework was done 
in Refs. lllOlJlin . including more exclusive sfudies of kinematical correlations. 

Very recently, the STAR collaboration measured for the first time transverse momen- 


Hli 


turn spectra of D* and mesons up to 6 GeV at = 200 11211 and 500 GeV m 
proton-proton collisions. Before, studies of charm production at RHIC through hadronic 
decay chaimels were performed only in Cu-Cu collisions lll4ll and in profon-profon scaf- 
fering but for D* mesons produced in jets 1151] . Up to now, the new STAR proton-proton 
data on open charm production were studied only in the context of high energy pA col¬ 
lisions in the Color Glass Condensate framework with the uninteerated eluon densities 

n 

from the solution of rcBK equation 11611 . 

Our aim here is fo make first theoretical analysis of the measured hadron-level charm 
differential cross section within the fcf-factorization approach. The open charm meson 
data allow us to make more direct comparison of the theoretical predictions and RHIC 
experimental results on heavy flavour production withouf including additional step re¬ 
lated to the semileptonic decays. Recently, the formalism of the U-factorization approach 
has been found fo give very good description of open charm IIitII and bottom ilSl. 11911 
production rafes and kinematical correlations in proton-proton scattering at ^/s = 7 TeV 
measured by the ALICE, ATLAS, CMS and LHCb experiments. However, a significant 
sensitivity of fheoretical predictions on fhe model of uninfegrafed (transverse momentum 
dependent) gluon distributions (UGDFs) used in calculations has been also reported. 

Therefore, if is very inferesting fo make similar study for the STAR experimental 
data on open heavy flavour production at = 200 and 500 GeV. This may be a 
good test of different models of UGDFs in the RHIC kinematical regime where one 
can probe parton (gluon) distributons at intermediate longitudinal momentum fractions 
x\lxj_ ^ 10“^ — 10“^. Here, we wish to pay particular attention on UGDF models 
favoured by the LHC data and on a new up-to-date parametrizations based on the HERA 
collider DIS high-precision data. The present study is an important extension of our pre¬ 
vious paper, where charm and bottom cross section at RHIC has been considered in the 
context of semilepfonic decays of open heavy mesons 


Ql. 
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Precise predictions for charmed mesons may also shed new light on the theoretical un¬ 
derstanding of non-photonic lepton production at RHIC. Our previous studies of these 
processes within the /cf-factorization approach were based on rather older models of 
UGDFs which may be the reason of the reported missing strength in description of the 
RHIC experimental data. Similar problem was also noticed within the NLO collinear 
calculations in the FONLL model, where only upper limits of the theoretical predictions 
are consistent with the relevant STAR and PHENIX data |^. Therefore, in the following 
paper we will also revise theorefical cross sections of the non-photonic lepton produc¬ 
tion at RHIC wihtin the fcrfactorizaton approach, taking as a reference point the results 
obtained in the analysis of the new hadron-level STAR data. In the following calcula¬ 
tion, except of updated UCDFs, we also take into account the effect of transformation of 
semileptonic decay functions between laboratory center-of-mass system) and rest 

frames of decaying D or B mesons. 


II. THEORETICAL FORMALISM 


Several different mechanisms play a role in heavy quark hadroproduction. In general, 
in the framework of QCD there are two types of the 0{a^) leading-order (TO) 2—7-2 
subprocesses: cjq ^ QQ and gg QQ 11200, often referred fo as heavy quark-antiquark 
pair creation. The first mechanism, ^^^-aImihilation, is important only near the thresh¬ 
old and at very large invariant masses of QQ system or extremely forward rapidities. 
This contribution is therefore especially important in the case of top quark production, 
however, for charm and bottom production, starting from RHIC, through Tevatron, up 
to the LHC it can be safely neglecfed. At high energies, production of charm and bot- 
fom flavours is dominated by the gluon-gluon fusion, which is the starting point of the 
following analysis. In the case of heavy quark production the 0(0.^) next-to-leading or¬ 
der (NLO) perturbative contributions have been found to be of special importance (see 
e.g. Ein i and have to be included in order to describe the heavy flavour high energy 
experimental data. 

The QED contributions, with one or two photons initiated reactions, such as ^g —?■ 
QQ> gl QQ, 77 —> QQ have been carefully studied in Ref. 11221] together with other 
sub-leading contributions to production of charm and were found fo be negligibly small 
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at high energies. 

In the following, the cross sections for charm and bottom quark production in proton- 
proton collisions are calculated in the framework of the /cf-factorization approach. This 
framework has been successfully applied for different high energy processes, including 


heavy quark production (see e.g. 1171] and references therein). According to this approach. 


the transverse momenta kt's (virtualities) of partons which initiate reaction are taken into 
account and the sum of transverse momenta of the final Q and Q no longer cancels. Then 
the TO differential cross section for the QQ pair production reads: 


da{hih2 QQX) 

dyidy2d^pi^td^P2,t 


E 

id 


d^ki^t d^k2^t 


1 


TZ TZ 16n^{XiX2sY 


\j^off-shel^2 
' g*g*-^QQ' 


( 2 . 1 ) 


X ^ [ki^t + ht - Pu - Pi,t) ^g{xi,k\,t) ^g{^2,'kit) / 


where Fg{xi, k\ and Tg{x 2 > ^ 21 ) UGDFs for the both colliding hadrons. The ex¬ 

tra integration is over transverse momenta of the initial partons. Explicit treatment of the 
transverse part of momenta makes the approach very efficient in studies of correlation 
observables. The two-dimensional Dirac delta function assures momentum conserva¬ 
tion. The unintegrated (transverse momentum dependent) gluon distributions must be 
evaluated at: 


^l,t ( \ I / \ 

Xx = ^exp(yi) + ^exp(y2), 

V S Y S 


X2 = ^exp(-yi) + ^exp(-y2), 
y S -y S 


where m/ t = y is the quark/antiquark transverse mass. In the case of heavy 

quark production at RHIC energies, especially in the central rapidity region, one test 
kinematical regime of x > 10“^. 

The TO matrix element squared for off-shell gluons is taken in the analytic form pro¬ 
posed by Catani, Ciafaloni and Hautmaim (CCH) 11231] . This analytic formula was basi¬ 
cally derived within the standard QCD framework and can be adopted to the numerical 
analyses. It was also checked that the CCH expression is consistent with those presented 
later in Refs. [l24]. 12511 arid in the limit oi k^^ —)■ 0, ^2 f Oh converges to the on-shell 
formula. 

The calculation of higher-order corrections in the fcf-factorization is much more com¬ 
plicated than in the case of collinear approximation. However, the common statement is 
that actually in the fcf-factorization approach with TO off-shell matrix elements some part 
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of higher-order corrections is effectively included. This is due to emission of extra gluons 
encoded in the unintegrated gluon densities. More details of the theoretical formalism 
adopted here can be found in Ref. ll^ . 

In the numerical calculation below we have applied several unintegrated gluon den¬ 
sities which are based on different theoretical assumptions. The Kimber-Martin-Ryskin 
(KMR) UGDF is derived from a modified DGLAP-BFKL evolution equation |26, 27] and 
has been found recently to work very well in the case of charm and bottom production 
at the LHG. A special emphasis here is put on the UGDF models obtained as a solution 
of GCFM evolution equation. Here we use an older Jung setBO parametrization II28I] and 
up-to-date JH2013 distributions ||29ll determined from the fits to HERA high-precision 
DIS measurements. The JH2013 setl is obtained from the fit to inclusive F 2 data only 
while JH2013 set2 is derived from the fit to both and F 2 data. The UGDFs based on 
BFKL and BK equations are not applied in the following analysis since they are dedicated 
to smaller-x values. 

In the calculation of charm and bottom quark cross sections the central value of nu¬ 
merical results is obtained with the renormalization and factorization scales }p- = }i\ = 

o Tfl^ I TFl^ 

Up = 2 and quark mass nic = 1.5 and = 4.75 GeV for charm and bottom, re¬ 

spectively. The uncertainties of the predictions are estimated by changing quark mass by 
±0.25 GeV and by varying scales by a factor 2. The gray shaded bands drawn in the fol¬ 
lowing figures represent these both sources of uncertainties summed in quadrature. The 
MSTW08LO 11301] collinear parton distribution function (PDF) is used for the calculation 
of the KMR unintegrated gluon density. 

The transition from quarks and gluons to hadrons, called hadronization or parton trag- 


models. In princi- 


mentation, can be so far approached only through phenomenologica 
pie, in the case of multi-particle final states the Lund string model 13IH and the cluster 
fragmentation model ]l32l] are often used. However, the hadronization of heavy quarks in 


non-Monte-Garlo calculations is usually done with the help of fragmentation functions 
(FF). The latter are similar objects as the parton distribution functions (PDFs) and provide 
the probability for finding a hadron produced from a high energy quark or gluon. 

Considering fragmentation of a high energy quark (parton) q with zero transverse 
momentum pt into a hadron q ^ h + X one usually assumes that the transition is soft 
and does not add any transverse momentum. In consequence it is a delta function in 
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transverse momentum. 

Defining D{z)dz as the probability for the quark q fragmenting into a hadron h which 
carries a fraction z of its longitudinal momentum, the spectrum of hadrons is given by 

ddfi 


dxhd'^pt,h 




dzdxD{z)-^S{x}i —zx) 

(A-X- 


or 


dcf}i 

dx}i 


*D(z)^ 

z dx 


( 2 . 2 ) 


(2.3) 


X=Xh/z 


This can be generalized to the fragmentation at finite pt with intrinsic transverse mo¬ 
mentum Kf. Thus the incoming quark has transverse momentum pt,q and the outgoing 
hadron transverse momentum can be decomposed as pt^h = K-t + zpt,q. Moving back to 
the situation where pt,q = 0 one can argue that the intrinsic transverse momentum Kt 
must be small. 

Introducing the transverse momentum dependent fragmentation probablity 
D(z, Kt)dzdKt one obtains the hadron spectum as 


do'll 


dxiid?-p 


t.h 


= / dzd^KtD{z,Kt)dXqd^pt,q- 


dCn 


-5{xii - ZXc,)P^\pt,h - k) - zpt,q 


'dXc,d^Pt,q 

However, one often neglects the intrinsic transverse momentum assuming: 

D{z,Kt) = D(z)P^^(Kt). 


Then the general formula reads 

do'll 


= / dzdxqd^pt,qD{z] 


do'n 


dxhd^Pih 
or integrating out the ^-functions 

do'n 

dxhd^Pih 


dxqd^pi^i 


-d{Xh - ZXq)P^\pt,li - Zpt, 


q) ' 


(2.4) 


(2.5) 


( 2 . 6 ) 


= IHd(z) 


do'n 


dXqd^Plq 


(2.7) 


Xq=Xj^/z 

Vt,q=Vt}id 


It is belived that this procedure provides correct implementation of small intrinsic trans¬ 
verse momentum into the splitting. Since the hadron on-shell four momentum is fully 
specified by Xh and pi^i, starting from the above formula one can calculate many different 
distributions. Especially conversion to the rapidity distributions is very simple because 
of the trivial jacobian: 

vdrr dn 

( 2 . 8 ) 


dxd'^pt dyd'^pi 
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It can be further written 


= fdzdxfy,,D(z)^^x,.5(x„-zx,)S^-^Kp,,,-zp,,) 


= I dz—(fpt,qD{z) ■ d{xh - zx^)d^^^ (pt^h - zpt,q). 


dXqd^Pt^q 


Neglecting masses, one has 


Xh = Xq = , 

VS VS 


( 2 . 10 ) 


so that. 


1^, = 1^. = 


= j dzdpqd^pt,qD{z) 5{yq - yh)d^^Hpt,h - zpt,q) 




( 2 . 11 ) 


vq=yh 

Vt,q=Vi,hl^ 


This way one gets (reproduces) the standard formula for the fragmentation in the case of 
light hadrons. 

When going to more general case, taking masses into account and introducing nii = 

\i P? + 


m? „ m? „ 

X], = Xq = , 

Vs 


Vh = log 




/ y<? = log 


Xq\[s 


( 2 . 12 ) 


(2.13) 


then, the z-dependent rapidity shift between quark and hadron reads 


= y<? - = log 


X\im[ 


= log 


(2.14) 


The delta function now becomes 


Xh&{X], - XqZ) = ^ 1^1 - ^ ^1 - j 

= ^ (l _ _ ^y _ /). 


( 2 . 15 ) 


9 



Finally, 


do'll 


dyhd^Pt,h 


= / dzdycjd^pt,cfD{z) 


d(Jn 


dyqd^Pt,q 


s(yh - y? + (Pi,h - zpt,q 


= ftD(z) 


dCn 


dycjd^p 


t,q 


( 2 . 16 ) 


yq=yh+^y 

Pt,q=Pt,h/^ 


In turn, using pt^q = Pt,h/^> rapidity shift Sy can be rewritten 


= 2 log 


pjh + K 

Vlh + 2^^? 


( 2 . 17 ) 


It is clear that the rescalling of the transverse momentum is the most important effect. 
This is because one deals with very steep functions of transverse momenta. From the 
reason that rapidity spectra are usually flat, or slowly varying, the shift 6y is not so im¬ 
portant. In fact, it is entirely negligible, if m^. The shift is most important at 

very small m^, nv^, where it becomes 


5y - log 


nih 

zma 




( 2 . 18 ) 


It is worth to notice, that at finite pt^h it should never be really large: small z is damped by 
the fact that the quark spectrum drops rapidly as a function of pt,q = Pt,h/‘Z- However, af 
Vt,h = 0/ that suppression causes an effect and the whole integral over z becomes impor¬ 
tant, with very small z causing large rapidity shifts. Fortunately, for heavy quarks, the 
fragmentation function is peaked at large z (see e.g. lISSll l. Moreover, one has to remem¬ 
ber, that taking into account the small-z region in the fragmentation function is theoreti¬ 
cally not warranted, since the standard DGLAP approach to fragmentation breaks down 
in this region. 

Taking all together, according to the above formalism, in the following numerical cal¬ 
culations the inclusive distributions of open charm and bottom hadrons h = D,B are 
obtained through a convolution of inclusive distributions of heavy quarks/antiquarks 
and Q ^ h fragmentation functions: 

da{pp —)• hhX) 


dyhd^Pt,}! 


‘ ^ QQx) 


/o z:- 


‘^yQd^Pt,Q 


( 2 . 19 ) 


yQ=yh 

Pt,Q=Pt,h/^ 


where pf q = ^ and z is the fraction of longitudinal momentum of heavy quark Q car¬ 
ried by a hadron h. The origin why the approximation typical for light hadrons assuming 
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that 1 /q is unchanged in the fragmentation process, i.e. = i/q, is also applied in the case 
of heavy hadrons was carefully clarified in the previous paragraph and is commonly ac¬ 
cepted. 

As a default set in all the following numerical calculations the standard Peterson 
model of fragmentation function 11341] with the parameters Sc = 0-02 and = 0.001 is 
applied. This choice of fragmenfation function and parameters is based on our previous 
theoretical studies of open charm production af fhe LHC illTI] . where defailed analysis 
of uncerfainties relafed wifh application of different models of FFs was done. Here, we 
decided nof fo repeaf all the previously analyzed issues and take into consideration only 
the most data-favoured scenario^. The main conclusions should not change when mov¬ 
ing from LHC fo RHIC energies and the uncertainties due to the fragmentation effects 
may be neglected with respect to those related to the perturbative part of the calculation. 

In the calculations of the cross sections for and D* mesons the fragmentation func¬ 
tions should be normalized to the relevant branching fractions BR(c —> D), e.g. from 
Ref. iSSt]. However, the measured by STAR differential distributions for D® and D* me¬ 
son are normalized to the parton-level cc cross section which simply means that the 
BR(c —^ D) = 1 should be taken in numerical calculations. 

Theoretical predictions for production of non-photonic leptons in proton-proton scat¬ 
tering is a three-step process. The whole procedure can be written in the following 
schematic way: 

dcr{pp —> l^X) dcr{pp —^ QX) 




( 2 . 20 ) 


dyed'^pt^e dyQd^Pt,Q 

where the symbol (8> denotes a generic convolution. Thus, the cross section for non- 
phofonic lepfons is a convolution of fhe cross section for heavy quarks wifh fragmenfa¬ 
tion function and with semileptonic decay function for heavy mesons. 

In principle, the semileptonic decays can be calculated 1361] . The simplest approach to 
describe the decays of D and B mesons is given by fhe specfafor model 1361] , where fhe 
QCD effecfs from the higher-order corrections between heavy Q and light q quarks are 
neglected. This model works better for bottom quarks since there the mass sufficiently 


^ This is also consistent with prescri^ion applied in the FONNL framework, where rather harder frag¬ 
mentation functions are suggested j9|]. 
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suppresses these corrections. In the case of charm, the QCD effects become more impor¬ 
tant but they can be also qualitatively modelled. 

Since there are many decay chaimels with different number of particles the above pro¬ 
cedure is not easy and rather labor-intensive. It introduces some model uncertainties and 
requires inclusion of all final state charmels explicitly. 


An alternative way to incorporate semileptonic decays in 


take relevant experimental input. For example, the CLEO BTH and BABAR 11381] collab¬ 


o theoretical model is to 


orations have measured very precisely the momentum spectrum of electrons/positrons 
coming from the decays of D and B mesons, respectively. This is done by producing res¬ 
onances: Y(3770) which decays into D and D mesons, and Y(4S) which decays into B 
and B mesons. 

This less ambitious but more pragmatic approach is based on purely empirical fits to 
(not absolutely normalized) CLEO and BABAR experimental data points. These electron 
decay functions should account for the proper branching fractions which are known ex¬ 
perimentally (see e.g. |j37H390). The branching fractions for various species of D mesons 
are different: 


BR(D+ ^ e+VeX) = 16.13 ± 0.20(stat.) ±0.33(syst.)%, 

BR(dO ^ e+VeX) = 6.46 ± 0.17(stat.) ± 0.13(syst.)%. (2.21) 

Because the shapes of positron spectra for both decays are identical within error bars we 
can take the average value of BR(D —> e Ve X) ~ 10% and simplify the calculation. In turn, 
the branching fraction of open bottom is found to be: 

BR(B ^eVeX) = 10.36 ± 0.06(stat.) ±0.23(syst.)%. (2.22) 

After renormalizing to experimental branching fractions the adjusted decay functions 
are then use to generate leptons in the rest frame of the decaying D and B mesons in a 
Monte Carlo approach. This way one can avoid all uncertainties associated with explicit 
calculations of semileptonic decays of mesons. 

In both cases the heavy mesons are almost at rest, so in practice one measures the 
meson rest frame distributions of electrons/positrons. With this assumption one can find 
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fit in lab frame (solid) boost in D rest frame (dashed) 




FIG. 1: Fits to the CLEO (left) and BABAR (right) data. The solid lines correspond to the 
parametrizations in the laboratory frames and the dashed lines to the meson rest frames, which 
represent incorporation of effects related to the non-zero motion of decaying mesons. 


a good fit to the CLEO and BABAR data with: 

/c?£o(p) = 12.55 (p + 0.02)2-55(0.98 - , (2.23) 

/b?bar(p) = (i 26.16 + 14293.09exp(-2.241og(2.51 - 0.97pf^ 

X (^-41.79 + 42.78exp(-0.5(|p - 1.27|)/1.8)®-28) . (2.24) 


In these purely empirical parametrizations p must be taken in GeV. 

In order to take into account the small effect of the non-zero motion of the D mesons 
in the case of the CLEO experiment and of the B mesons in the case of the BABAR experi¬ 
ment, the above parametrizations of the fits in the laboratory frames have to be modified. 
The improvement can be achieved by including the boost of the new modified rest frame 
functions to the CLEO and BABAR laboratory frames. The quality of fits from Eqs. (12.231) 
and (|2.24|) will be reproduced. The D and B rest frame decay functions take the following 
form: 


ftiioip) = 12.7(p + 0.047)2-22(0.9 - p)2-2i , (2.25) 

fBABARip) = (l26.16 +14511.2 exp(-1.931og(2.7-1.0825p)2j 

X (^-41.79 + 42.78exp(-0.5(|p - 1.27|)/1.8)®-28) . (2.26) 

Both, laboratory and rest frame parametrizations of the semileptonic decay func¬ 
tions for D and B mesons are drawn in Eig. [T] together with the CLEO (left panel) and 
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BABAR (right panel) experimental data. Some small differences between the different 
parametrizations appear only at larger values of electron momentum. The influence of 
this effect on differential cross sections of non-photonic leptons is expected to be negli¬ 
gible and will be shown when presenting numerical results. Our analytical formulas for 
the rest frame decay functions only slightly differ from those obtained in Ref. 11401] . 

The theoretical model for non-photonic lepton production in hadronic reactions de¬ 
scribed here has been recently found to give a very good description of the experimental 
data collected with the ALICE detector at the LHC |41,42]. 


III. NUMERICAL RESULTS 

The total cross section for charm production in pp scattering extracted from the STAR 
measurement of and D* mesons at ^/s = 200 GeV is dec = 797 ± 2 IO 295 Corre¬ 
sponding calculated total cross sections is = 541 with the CCFM Jung setBO UGDR 
The calculated value is consistent with the measured value taking into account large ex¬ 
perimental uncertainties. The total cross section at ^/s = 500 GeV using the same UGDF 
is predicted to be tr®!™ = 1006 pb. 

The STAR collaboration also carried out a measurement of charm production cross 
sections at midrapidity Gomparison of the experimental results and the theo¬ 

retical ones is presented in Table HI Here again the GGFM Jung setBO UGDF give results 
consistent with the measurements. 

TABLE I: The midrapidity cross section for charm production in proton-proton scattering 

at ^/s = 200 and 500 GeV: the STAR results versus results of calculation with the Jung setBO 


UGDF. 


Experiment: STAR, ^|y=o 

Theory: Jung setBO 

x/s = 200 GeV 

170 ± pb 

130 pb 

x/s = 500 GeV 

217 ± 86 ± 73 pb 

191 pb 
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A. Open charm mesons 


Figure |2] presents transverse momentum distributions of charmed mesons in proton- 
proton collisions at ^/s = 200 GeV for |i/d| < 0.5 (left panel) and at = 500 GeV for 
|i/d| < 1-0 (right panel) together with the STAR data points. Both experimental data 
and theoretical results are normalized to the cc parton-level cross section dividing by the 
c —?■ D*, fragmentation fractions. Results of numerical calculations obtained with the 
KMR (dotted line), the JH2013 setl (long-dashed-dotted line), set2 (long-dashed line) and 
the Jung setBO (solid line) UGDFs are shown. At both energies very good description of 
the experimental data is obtained with the Jung setBO UGDF Results calculated with the 
JH2013 setl UGDF overestimate the data points in the whole range of measured pt's. The 
JH2013 set2 and the KMR UGDFs significantly underestimate the distribution measured 
at = 200 GeV and the situation is only slightly improved in the case of the higher 
energy where the results of both of them reach the two last data points at larger transverse 
momenta. 




FIG. 2: The transverse momentum distribution of D mesons normalized to the parton-level cross 
section at ^/s = 200 (left) and 500 GeV (right). The STAR experimental data points are compared 
to the results of the fcf-factorization calculations with the KMR (dotted line), the JF12013 setl (long- 
dashed-dotted line), set2 (long-dashed line) and the Jung setBO (solid line) UGDFs. 

Main uncertainties of the theoretical calculations coming from the perturbative part 
are shown in Fig.|3l The shaded bands represent uncertainties of the calculations with the 
Jung setBO UGDF related to the choice of the factorization and/or renormalization scales 
and those due to the charm quark mass. The result from the FONLL approach is also 
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drawn for comparison. The uncertainties are larger at lower transverse momenta, where 
the effects of quark mass uncertainties are more important, and decrease with increasing 
pt- The FONLL predictions underestimate the experimental data almost in the whole 
measured range. Their central value coincides with the lower-limit of the fcf-factorization 
predictions with the Jung setBO UGDR 




FIG. 3: The uncertainties of the theoretical predictions for Jung setBO UGDF at y/s = 200 (left) and 
500 GeV (right). Uncertainties due to the choice of the factorization and/or renormalization scales 
and those related to the charm quarks mass are summed in quadrature. In the case of ^/s = 200 
GeV data the results obtained within the FONLL framework are drawn for comparison. Details 
are specified in the plots. 


B. Non-photonic electrons 


A first theoretical investigation of the non-photonic electron production at RHIC 
within the framework of the /cf-tactorization was performed in Ref. Ilj. Some missing 
strength in the description of the measured differential distributions has been reported 
there, especially in the region of small transverse momenta. In the meantime, the STAR 
collaboration has published new measurements of non-photonic electrons with separated 
charm and bottom contributions [71]. Therefore, it is very interesting to make a revision 
of theoretical cross sections taking into account the new results from the hadron-level 

n 

analysis of charm production at RHIC. Here, our previous results from Ref. piflll are up¬ 
dated by the application of the Jung setBO UGDF which as was shown in the previous 
subsection works very well for the STAR data on open charm meson production. 
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The experimental cross sections for charm and bottom flavoured electrons measured 
at RHIC are collected in Table im The values calculated with the Jung setBO UGDF are 
consistent with the measurements. 

TABLE II: The experimental and theoretical cross sections for non-photonic electron pro¬ 


duction in proton-proton scattering at \/s = 200 GeV. 


Experiment 

Theory, Jung setBO 

STAR, pp —)■ ccX —)■ eX' 

6.2±0.7±1.5nb 

g:li,,.o,3<px <10GeV 

7.55 nb 

STAR, pp bbX eX' 

4.0±0.5±l.lnb 

^U=o,3<p^<10GeV 

6.65 nb 

PHENIX, pp ccX eX' 

5.95 ± 0.59 ± 2.0 ph 

^U=o,p^>0.4GeV 

5.09 ph 


Figure m shows the transverse momentum distributions of electrons from semileptonic 
decays of charm flavoured hadrons He (left panel) and from bottom hadrons Hi, (right 
panel) measured by STAR. The experimental data is compared to the numerical results 
for the Jung setBO UGDF. The theoretical uncertainties coming from the perturbative part 
of calculations are also shown for completeness. The rest frame semileptonic decay func¬ 
tions from Eqs. (I2.25D and (I2.26D are used. It is also assumed that the charm and bottom 
baryons decay semileptonically in the same way as D and B mesons, and therefore bary- 
onic contributions may be effectively included by treating the baryons as mesons and 
taking BR(c —> D; & —^ B) = 1. The numerical results very well descibe the experimental 
data. The central value of the Jung setBO UGDF give distributions that are sligthly above 
the predictions of the FONLL central value, especially in the small-pt region. In this case 
also the JH2013 setl UGDF reasonably describes the data points taking into account ex¬ 
perimental uncertainties. As in the case of open charm data, the lines that corespond to 
the KMR and JH2013 set2 UGDFs lie much below the measured lepton distributions for 
both, charm and bottom components. 

In Fig. 13 the results for summed contributions of charm and bottom flavours are 
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FIG. 4: Transverse momentum distributions for electrons from semileptonic decays of charm (left) 
and bottom hadrons (right) measured in pp scattering at = 200 GeV. The STAR experimental 
data are compared to the fcf-factorization theoretical predictions obtained with different UGDFs 
as well as to the FONLL results. Theoretical uncertainties due to quark mass and scales variation 
are also shown. Further details are specified in the figures. 




FIG. 5: Transverse momentum distributions of electrons coming from both charm and bottom 
hadrons summed together Hc+f,, measured in pp-scattering at \/s = 200 GeV. The STAR and 
PFIENIX experimental data are compared to the theoretical predictions obtained with the Jung 
setBO UGDF. Separated charm and bottom contributions (left) and theoretical uncertainties due 
to quark mass and scales variation (right) are also shown. The FONLL predictions are drawn for 
comparison. Lurther details are specified in the figures. 


shown. Here, the results of calculations are compared to the experimental distributions 
of heavy flavour electrons, that contain both charm and bottom components. The left 
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panel presents separately charm and bottom contributions as well as the sum of them. In 
the right panel, the uncertainties of the predictions for the Jung setBO UGDF are drawn 
together with the lines corresponding to the FONLL results. In contrast to the previous 
studies in Ref. 11101] , here the fcf-factorization results give excellent description of the STAR 
and PHENIX data, sligthly better than those from the FONLL approach, which are almost 
identical to the lower for the Jung setBO UGDF. The crossing point between charm and 
bottom components is found fo lies roughly af pt = 4 GeV, which is in agreemenf with 
other theoretical investigations (see e.g. Ref. fli). 




FIG. 6: The effect of using two different sets of the semileptonic decays functions, discussed in 
the present paper, at -\/s = 200 GeV. The results obtained with the laboratory frame fits from 
Eqs. (|2.23D and (12.241) are compared to the meson rest frame fits from Eqs. (12.251) and (|2.26|) . 


Finally, fhe effecfs relafed fo fhe fitting procedure of the GLEO and BABAR semilep¬ 
tonic data are depicted in Fig. [6l Here, both sets of the semileptonic decay functions are 
used. As can be observed from the figure fhese effecfs do nof really affecf the electron 
spectra at RHIG. The difference is very small and sligthtly increases at higher transverse 
momenta. In the case of charm flavour (leff panel), the application of the boosted decay 
function leads to a damping of fhe cross section by abouf 20% af pt = 10 GeV. For fhe 
bottom flavour (right panel) the corresponding suppresion is only about 5%, which is 
completely negligible. 
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IV. CONCLUSIONS 


In this paper we have dicussed production of charm mesons and non-photonic elec¬ 
trons in proton-proton scattering at the BNL RHIC. The calculation of the charm quark- 
antiquark pairs has been performed in the framework of fcf-factorization approach which 
effectively includes higher-order pQCD corrections. 

We have used different models of unintegrated gluon distributions from the literature, 
including those that were applied recently to describe charm data at the LHC and others 
used to descibe HERA deep-inelastic scattering data. The hadronization of heavy quarks 
to mesons has been done by means of fragmentation function technique. The theoretical 
transverse momentum distributions of charmed mesons has been compared with recent 
experimental data of the STAR collaboration collected at ^/s = 200 and 500 GeV. We have 
carefully quantified uncertainties related to the choice of factorization/renormalization 
scales as well as quark/antiquark masses. We have obtained very good agreement with 
the measured cross sections for the Jung setBO UGDE Furthermore, our results have been 
compared with the results of the FONLL model. The two approaches give rather similar 
results. 

Semileptonic decays of charmed and bottom mesons have been included via empirical 
decay functions fitted to the GLEO and BABAR {e^e~) data for vector meson decays. We 
have shown that the inclusion of kinematical boost from meson (D or B) rest frame to the 
e~^e~ center of mass (laboratory) system leads to only small modifications of the resulting 
decay functions and as a consequence also for the distributions of non-photonic elec¬ 
trons in proton-proton collisions at the RHIG energies. Gonsequently we have obtained 
a rather good description of the electron/positron transverse momentum distributions 
of the STAR collaboration with the same UGDF as for the charmed mesons. This also 
demonstrates indirectly consistency of the meson and non-photonic electron data. 
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